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Abstract. - We present a detailed bifurcation structure and associated flow patterns for low- 
Prandtl number (P = 0.0002, 0.002, 0.005, 0.02) Rayleigh-Benard convection near its onset. We 
use both direct numerical simulations and a 30-mode low-dimensional model for this study. We 
observe that low-Prandtl number (low-P) convection exhibits similar patterns and chaos as zero- 
P convection [1], namely squares, asymmetric squares, oscillating asymmetric squares, relaxation 
oscillations, and chaos. At the onset of convection, low-P convective flows have stationary 2D rolls 
and associated stationary and oscillatory asymmetric squares in contrast to zero-P convection 
where chaos appears at the onset itself. The range of Rayleigh number for which stationary 
2D rolls exist decreases rapidly with decreasing Prandtl number. Our results are in qualitative 
agreement with results reported earlier. 
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Thermal convection is an important problem of classical 
physics, engineering, geophysics, atmospheric physics, and 
astrophysics. As a consequence of its inherent nonlinear- 
ity, good understanding of convection eludes us even after 
sustained efforts for more than a century. Convection in an 
arbitrary geometry is quite complex, so researchers have 
focusscd on convection between two conducting parallel 
plates. This idealized problem is called Rayleigh-Benard 
convection (RBC). 

Convective flow in RBC is characterized by two non- 
dimensional parameters: the Rayleigh number R which is 
the ratio of the destabilizing force to the stabilizing force, 
and the Prandtl number P which is the ratio of the kine- 
matic viscosity to the thermal diffusivity. RBC shows a 
wide range of phenomena including instabilities, patterns, 
chaos, spatio-temporal chaos, and turbulence for various 
ranges of R and P [2]. For low-Prandtl number (low-P) 
and zero-Prandtl number (zero-P) convection, the iner- 
tial term u • Vu becomes quite important and generates 
vertical vorticity. As a result, the flow pattern becomes 
three-dimensional, and oscillatory waves along the hori- 
zontal axes are generated just near the onset of convec- 
tion. On the contrary, for large-Prandtl number (large-P) 
convection, vertical vorticity is absent near the onset and 
the two-dimensional (2D) rolls survive till large Rayleigh 



numbers [3]. 

In this letter we investigate the bifurcation scenario for 
low-P convection and compare it with the recent zero-P 
convection results of Pal et al. [1]. We will also explore the 
origin of chaos in low-P convection. The relevant equa- 
tions (under Boussincsq approximation) in nondimension- 
alizcd form arc 



dt(V 2 v 3 ) 
d t uj 3 

P(d t 6 + (v ■ V)0) 
V • v = 



V 4 w 3 + RV 



e 3 • V x [(u • V)v - (v • V)i4l) 
[(w • V)w 3 - (v • V)w 3 E2) 
(3) 
(4) 



V 2 w 3 



v 3 + w 2 e, 



where v = (v\, 1)2, V3) is the velocity field, 9 is the devia- 
tion in the temperature field from the steady conduction 
profile, u> = V x v is the vorticity field, R — ag(AT)d 3 /isk 
is the Rayleigh number, P = v/k is the Prandtl num- 
ber, e»3 is the unit vector along the vertical (buoyancy) 
direction, and V 2 ^ = d xx + d yy is the horizontal Lapla- 
cian. Here, v and k arc fluid's kinematic viscosity and 
thermal diffusivity respectively, d is the vertical distance 
between the two plates, and AT is the temperature differ- 
ence between the horizontal plates. For the above nondi- 
mensionalization, we used d as the length scale, d 2 jv as 
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the time scale, and vAT/k as the temperature scale. We 
consider perfectly conducting and free-slip conditions at 
the top and bottom plates and periodic boundary condi- 
tions along the horizontal directions. Wc will mostly use 
the reduced Rayleigh number r = R/ R c , where R c is the 
critical Rayleigh number, as the control parameter for our 
bifurcation analysis. 

Instabilities and chaos in low-P and zero-P convection 
have been widely studied using analytical tools (e.g., ap- 
proximation techniques), experiments, and numerical sim- 
ulations. Busse [4] showed using perturbative analysis 
that for small Prandtl numbers, the 2D rolls become un- 
stable as a result of oscillatory three-dimensional distur- 
bances when the amplitude of the convective motion ex- 
ceeds a critical value. Busse and Bolton [5] argued that 
under free-slip boundary conditions, the stability of the 2D 
rolls exists only for Prandtl numbers above a critical value 
P c , which is around 0.543. Clever and Busse [6] extended 
the oscillatory instability analysis to no-slip boundary con- 
ditions and showed that the convective rolls are unstable 
for Prandtl numbers less than about 5. 

Siggia and Zippelius [7] analyzed convective instabili- 
ties using amplitude equations and reasoned that for finite 
Prandtl numbers, the vertical vorticity modes are impor- 
tant and the flow is three-dimensional. Fauve et al. [8] 
investigated the origin of instabilities in low-P convection 
using the phase dynamical equations and argued that the 
instability always saturates into travelling waves as ob- 
served in experiments and numerical simulations. Jenkins 
and Proctor [9] studied the transition from 2D rolls to 
square patterns in RBC for different Prandtl numbers us- 
ing analytical tools. Jones et al. [10] studied connection 
of the Prandtl number and the Nusselt number using per- 
turbative techniques and numerical simulations. Kumar 
et al. [11] investigated zero-P convection using a 6-mode 
model and showed the existence of wavy instability at the 
onset. They demonstrated that the saturation of the grow- 
ing nonlinear roll solutions is due to the wavy instability. 

Interesting experiments have been performed to explore 
the instabilities and chaos in convective flows near the on- 
set. Willis and Deardorff [12] carried out experiments on 
air (P = 0.7) and observed time dependence as a result 
of wavy instability. Rossby [13] demonstrated the exis- 
tence of wavy instability in mercury (P = 0.02). Krish- 
namurti [14] performed extensive convection experiments 
on a variety of fluids (P = 0.1 — 50) and observed station- 
ary and time-dependent three-dimensional patterns in her 
experiments; time dependent rolls appear at much lower 
Rayleigh numbers for low-P convection compared to that 
for large-P convection. 

Numerical simulations complement experiments in the 
investigation of instabilities and chaos in convection. 
Lipps [15] simulated convective flows in air (P w 0.7) 
and investigated various states including 2D rolls, time- 
periodic and aperiodic convection as a function of Rayleigh 
number and aspect ratio. Bolton and Busse [16] studied 
the stability of the steady convection rolls with respect 



to arbitrary three-dimensional infinitesimal disturbances 
and observed stability of steady solutions for only a small 
regime of Rayleigh numbers and wavenumbers. Clever 
and Busse [17] simulated convective flows for P = 0.02 
with rigid boundary conditions using the Galerkin tech- 
nique and observed asymmetric travelling wave patterns. 
Mcncguzzi et al. [18] performed three-dimensional convec- 
tive simulations for P = 0.2 fluid under stress-free bound- 
ary conditions, and for P = 0.025 for no-slip boundary 
conditions. For P = 0.2 they observed stationary, peri- 
odic, biperiodic and chaotic regimes as the Rayleigh num- 
ber is increased. For the latter, they found stationary 
and time-periodic solutions only. Thual [19] studied zero- 
P and low-P convection using numerical simulations for 
both no-slip and free-slip boundary conditions. In his sim- 
ulations he observed supercritical oscillatory instabilities, 
competition between two-dimensional rolls, square and 
hexagonal patterns, travelling and standing waves, and 
chaos for both zero-P and low-P convective flows. Ozoe 
and Hara [20] performed numerical simulations for Prandtl 
numbers P = 0.001 — 0.1 in a two-dimensional rectangular 
enclosure of aspect ratio 4 and observed that the time- 
dependent behaviour sets in at the onset of convection for 
P = 0.01, thus decreasing the value of P c to much lower 
limit than that predicted by Busse and Bolton [5] . 

Recently Pal et al. [1] performed numerical simulations 
and bifurcation analysis using a low-dimensional model 
of zero-P convection under free-slip boundary conditions. 
They could explain the origin of squares (SQ), asymmetric 
squares (ASQ), oscillating asymmetric squares (OASQ), 
relaxation oscillations with an intermediate square pat- 
tern (SQOR), and chaos using a series of bifurcations. 
Pal et al. [1] start their analysis from square patterns 
at r = 1.4. Under reduction of r, SQ bifurcates into 
asymmetric squares through a pitchfork bifurcation, which 
in turn transforms into oscillating asymmetric squares 
through a Hopf bifurcation. OASQ transforms to SQOR 
that subsequently becomes chaotic. The chaotic attractor 
continues till the onset of convection. There have been 
other attempts to study RBC using low-dimensional mod- 
els for other ranges of Prandtl numbers [21]. 

In this letter we perform numerical simulations and bi- 
furcation analysis of low-P convection along similar lines 
as Pal et al. [1] to investigate the origin of patterns and 
chaos. We also carry out a comparative study of low- 
P convection with zero-P convection and show how the 
limiting behavior of zero-P convection is obtained as P 
approaches zero. 

Wc perform direct numerical simulations (DNS) of con- 
vective flows for P = 0.02,0.005,0.002 and 0.0002 for r 
ranging from 1 to 1.247 on 64 3 rectangular grid. We use a 
pseudo-spectral procedure with fourth-order Runge-Kutta 
as the time-stepping scheme. The aspect ratio of our sim- 
ulations is 2^2 : 2 V2 : 1. Eqs. (Q]) and © provide us 
an estimate of dt for the DNS (dt ~ P6/v 3 ~ P/R). For 
low-P simulations this dt is typically much smaller than 
that obtained using the CFL condition. For example, for 



Patterns and bifurcations in low-Prandtl number Rayleigh-Bcnard convection 



P = 0.0002, dt ~ that makes numerical simulations 
very demanding. In DNS we observe stationary 2D rolls, 
stationary asymmetric square patterns (ASQ), wavy rolls, 
oscillatory asymmetric squares (OASQ), relaxation oscil- 
lation with squares (SQOR), and symmetric square pat- 
terns (SQ) for P = 0.02. A detailed description of the 
above phenomena in terms of bifurcation diagrams is de- 
scribed below. 

For a detailed bifurcation analysis of this problem, we 
construct a 30-mode low-dimensional model using the en- 
ergetic modes of the DNS [1]. We pick 11 large-scale verti- 
cal velocity modes: Wioi,Woii, W112, W2111 W121, W301, 
W103, W013, W031, W202, W022; 12 large-scale modes: 

0101! 0011 i 0112i 0211, 0121, 0301, 0103, 0013, 0031, 0202, 0O22i 

6*002; and 7 large-scale vertical vorticity modes: Zuq, ^112, 
■^211, ^121 1 ^220! ^310, ^130 • The three subscripts here are 
the indices of the wavenumbers along the x, y, and z di- 
rections. All the Fourier modes, 30 in total, are taken 
to be real. This model is inspired by the 13-mode model 
for zero-P convection [1] . Next we investigate the origin of 
the convective flow patterns using a bifurcation diagram of 
this low-dimensional model following the approach of Pal 
et al. [1]. The fixed point(s) are generated numerically us- 
ing the Newton- Raphson method for a given r, and they 
are subsequently continued using a fixed arc-length based 
continuation scheme [22] for the neighbouring r values. 
The stability of the fixed points and the parameter val- 
ues corresponding to bifurcations are ascertained through 
an eigenvalue analysis of the Jacobian. New branches of 
steady solutions (fixed points at pitchfork, while periodic 
solutions at Hopf ) are generated by calculating and contin- 
uing these solutions. The stability of these new solutions 
is ascertained analogously to obtain subsequent bifurca- 
tions and new branches of solutions. A MATLAB package 
named MATCONT was used for some of the above anal- 
ysis. 

The state space or the phase space of our low- 
dimensional model is 30 dimensional. However among all 
the modes, W101, Won, 0ioi, an d 0on are the most impor- 
tant ones; the behaviour of the rolls along the x and the 
y directions depends critically on these modes. In Fig. Q] 
we plot the numerical values of the fixed points W101 and 
Won for P = 0.02 as a function of r. The stable fixed 
points are shown as solid lines, while the unstable ones 
are shown as dotted and dashed lines. At r = 1, the 
conduction state (cyan curve) becomes unstable and bi- 
furcates into four stable 2D rolls (shown as purple curves) 
and four unstable symmetric squares (SQ, shown as black 
dotted lines) through a codimension-2 pitchfork bifurca- 
tion associated with double zero eigenvalues. The stable 
2D rolls along the x axis have Won = 0, while the rolls 
along the y axis have W101 = 0. The unstable symmet- 
ric square (SQ) solutions satisfy |Wioi| = |Woii|. With 
an increase of r, the stable 2D rolls bifurcate into stable 
patterns called asymmetric squares (ASQ) shown as the 
solid blue lines in Fig. [TJ These ASQ patterns lose sta- 
bility at some point through a Hopf bifurcation. With 



a further increase in r, the unstable branch of ASQ so- 
lutions (dotted blue curves) regains stability through an 
inverse Hopf bifurcation resulting in stable ASQ solutions 
(solid blue curves). These stable ASQ solutions subse- 
quently meet the unstable SQ patterns (represented by 
dotted black curves originating from r = 1) in yet an- 
other inverse bifurcation, and stable SQ patterns (solid 
black curves) are formed. These stable SQ patterns sub- 
sequently lose stability beyond r — 1.247 and complex 
chaotic attractors are generated, however an analysis in 
this regime is beyond the scope of this letter. This se- 
quence of fixed points remains the same for all finite but 
small values of P. However, the range of r corresponding 
to these solutions changes with P. 
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Fig. 1: For P = 0.02, the three dimensional view of the bi- 
furcation diagram computed from the low-dimensional model 
exhibiting stable fixed points (solid curves), and unstable fixed 
points (dotted and dashed curves). The purple solid and 
dashed lines represent the stable and unstable 2D roll solutions 
respectively. The black, blue, and cyan curves depict station- 
ary squares (SQ), asymmetric stationary squares (ASQ), and 
conduction state respectively. 

We will describe the above mentioned fixed points and 
associated time-dependent patterns using a bifurcation di- 
agram in the range of 0.95 < r < 1.247. Figure [2] illus- 
trates the bifurcation diagram for P = 0.02 where we plot 
the positive value of (W101 )extremum as a function of r. At 
r = 1 the conduction state bifurcates to stationary 2D rolls 
(purple curve) through a codimcnsion-2 pitchfork bifurca- 
tion. As described above unstable stationary SQ (shown 
as dashed black curve) with W101 = ±Won are also gen- 
erated as a result of this bifurcation. As r is increased 
further, at r k 1.0035, the stationary 2D rolls bifurcate 
to ASQ patterns through another pitchfork bifurcation 
(solid blue curves). Subsequently, at r sa 1.0114, ASQ 
patterns bifurcate to limit cycles (red curves) through a 
Hopf bifurcation. These limit cycles represent oscillatory 
asymmetric square patterns (OASQ). Note that the limit 
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observed near the onset of convection. The chaotic attrac- 
tors Chi, Ch2, and Ch3 of zero-P convection [1] are absent 
for P = 0.02. In the following discussion we will con- 
sider the bifurcation diagrams for lower Prandtl numbers 
(P = 0.005,0.002,0.0002) for which chaotic attractors are 
observed. 

We performed DNS for P = 0.02 and observed similar 
behaviour as the low-dimensional model. A close similar- 
ity between the model and the DNS is due to the fact 
that the modes for the low-dimensional models were cho- 
sen from DNS runs by identifying the energetic modes. 
This exercise is possible near the onset of convection where 
only a limited number of modes are excited. Comprehen- 
sive analysis of DNS and low-dimensional model will be 
presented elsewhere. 



Fig. 2: Bifurcation diagram of the low-dimensional model for 
P = 0.02 in the range 0.95 < r < 1.247. The stable branches 
corresponding to conduction state, 2D rolls, stationary squares 
(SQ) and stationary asymmetric squares (ASQ) are represented 
by solid cyan, solid purple, solid black, and solid blue curves 
respectively. The red and green curves depict the extrema of 
oscillatory asymmetric squares (OASQ) and relaxation oscil- 
lation with intermediate square patterns (SQOR) respectively. 
The dashed curves represent unstable solutions. 

cycles are quite close to the axes. The oscillatory pat- 
terns corresponding to these limit cycles form standing 
waves along the roll axis that have been discussed earlier 
by Thual [19]. The time-period and amplitude of these 
limit cycles increase as r is increased. At r k 1.0184, 
these limit cycles form homoclinic orbits of the unstable 
SQ fixed points (saddles) originating from r = 1. 

Beyond r ~ 1.0184, there is a smooth transition from 
homoclinic orbits to regular limit cycles that correspond to 
relaxation oscillations between two roll solutions along the 
x and the y directions with intermediate square patterns 
(SQOR, illustrated as green curves in the bifurcation di- 
agram). The SQOR patterns and all subsequent patterns 
generated at higher r values are similar to that observed 
by Pal et al. [1] for zero-P convection. In brief, the limit 
cycles corresponding to SQOR form another set of homo- 
clinic orbits at r ss 1.1034 after which each of the homo- 
clinic orbits bifurcate into two separate limit cycles. These 
limit cycles (OASQ shown as red curves) diminish in size 
as r is increased until they transform to stable fixed points 
(ASQ shown as blue curves) through an inverse Hopf bi- 
furcation. As r is increased further, at r ~ 1.183, the 
stable ASQ branches meet the unstable SQ branch (the 
dashed black curve) resulting in stabilization of the SQ 
pattern (solid black curves) as a consequence of an inverse 
pitchfork bifurcation. Thus, convective flow for P = 0.02 
exhibits SQ, ASQ, OASQ, and SQOR patterns that are 
common to zero-P convection [1]. In addition, stationary 
2D rolls and associated stationary and oscillating ASQ are 
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Fig. 3: Bifurcation diagram of the low-dimensional model for 
P = 0.005 in the range 0.95 <r< 1.247. The color convention 
is same as that for P = 0.02 (Fig. [2j. The chaotic attractors 
are shown in brown colour. A zoomed view of the bifurcation 
diagram for the chaotic regime is shown in the inset where the 
x axis is chosen as log(r — 1) to highlight the behaviour near 
r = 1. The three attractors Chi, Ch2 and Ch3 are shown as 
(i), (ii), and (iii) respectively. 

After the discussion on the bifurcation scenario for P = 
0.02, we turn to lower Prandtl numbers P = 0.005, 0.002, 
and 0.0002. The bifurcation diagram for P = 0.005, shown 
in Fig. [3l is similar to the corresponding figure for P = 0.02 
with a major difference that for P = 0.005, the chaotic at- 
tractors Chi, Ch2, and ChS of zero-P convection [1] appear 
near the onset in the band of r = 1.000685—1.0068. As we 
increase r from 1, we observe 2D straight rolls, ASQ, and 
OASQ just like P = 0.02, however, at lower r values (see 
Tabled]). The phase space projection on the Wioi-Won 
plane of two of the limit cycles corresponding to OASQ 
are shown in figure [Ha). Subsequently these limit cycles 
appear to approach their basin boundary (the horizontal 
axis of the figure), and the system becomes chaotic. This 
phenomenon is due to a "gluing bifurcation" [23], and it 
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could be related to the "attractor merging crisis" [24] . The 
resulting chaotic attractor, whose phase space projection 
is shown in Fig.[^b), is same as the Ch3 attractor of zero- 
P convection [1]. As r is increased further, the various Ch3 
attractors merge together in yet another crisis to form a 
single large chaotic attractor Ch2. At a larger r value, the 
Ch2 attractor breaks into four separate chaotic attractors 
Chi. Phase space projections of the chaotic attractors Ch2 
and Chi for P = 0.005 are shown in Figs. IDJc) and IU(d) 
respectively. The Chi chaotic attractors become regular 
for r > 1.0068 giving rise to the SQOR limit cycles. The 
subsequent patterns and the bifurcation diagram are same 
as those for P = 0.02. Note that the Chi chaotic attrac- 
tors are generated as a result of "homoclinic chaos" [1]. 



Table 1: Table depicting the reduced Rayleigh number r p at 
which the 2D rolls bifurcate to ASQ, r t where ASQ patterns 
bifurcate to limit cycles, and the frequency lu of the limit cycle 
at the Hopf bifurcation point. These values are computed for 
Prandtl numbers P = 0.02, 0.005, 0.002 and 0.0002. 
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The bifurcation scenario presented above for P = 0.005 
has all the features (including chaos) of zero-P convection 
along with the stable 2D rolls and associated ASQ and 
OASQ near the onset. Note that for P = 0.02, the window 
of chaos ( Chi to Ch3) is absent as the homoclinic intersec- 
tion point of the left OASQ limit cycle precedes the gluing 
bifurcation point. The patterns and chaotic attractors ob- 
served in the low-dimensional model for P = 0.005 are also 
observed in DNS. Phase space projections of the chaotic 
attractors ( Ch3 and Chi) obtained in DNS arc shown in 
Figs. (Ha) (for r = 1.0023) andGJb) (for r = 1.0068) re- 
spectively. The chaotic attractors obtained from the DNS 
data are not thick, however an enlarged view of the at- 
tractor in the inset shows that the attractor is chaotic. 

When we lower the value of the Prandtl number even 
further from P = 0.005, the stationary 2D rolls and the 
associated ASQ patterns occur even nearer to r = 1 as 
evident from the entries of Table [T] It can be noted from 
FigUthat the windows of 2D rolls and ASQ for P = 0.0002 
is very small. The inset of Fig. [5] shows that the 2D 
rolls and ASQ occur for 1 < r < 1 + 2.5 x 10 -7 and 
1 + 2.5 x 10~ 7 < r < 1 + 1.8 x 10 -7 respectively. Con- 
sequently it would be very difficult to observe these 2D 
rolls and the related ASQ patterns in experiments and 
in DNS; a small experimental or numerical noise will be 
sufficient to push the system out of this narrow region to 
the chaotic region Chl-Ch3, which is a feature of zero-P 
convection. Thus our study indicates that zero-P convec- 
tion is an appropriate limiting case of low-P convection in 
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Fig. 4: Phase space projections of the attractors on the Wioi- 
ITon plane obtained from the low-dimensional model for P — 
0.005 near r = 1: (a) two of the limit cycles at r = 1.0006844; 
(b) the chaotic attractor Ch3 at r = 1.0023; (c) the chaotic 
attractor Ch2 at r = 1.0053; (d) the chaotic attractor Chi at 
r = 1.0064, with the inset depicting the chaotic nature of the 
attractor. 

terms of bifurcations near the onset. For practical pur- 
poses, the bifurcation diagram for very low-P convection 
is topologically equivalent to zero-P convection. 

It can be seen from Table [1] that the branch points cor- 
responding to ASQ (r p ) and OASQ (r t ) asymptotically 
approach r = 1 as P — > 0. In Table [T] we also list the 
imaginary part of the eigenvalue (w) of the stability ma- 
trix at the Hopf bifurcation (r = rt). We observe that 
r p — 1, rt — 1, and lu appear to vary approximately as P 2 
with prefactors around 10, 30 and 3 respectively. The P 2 
dependence is consistent with the theoretical predictions 
of Busse [4] that r t — 1 > 0.310P 2 for free-slip boundary 
conditions. However the multiplying factor of Busse and 
our model differ by an order of magnitude. Also, Busse 
and Bolton [5] predict an absence of stable 2D rolls for 
P < P c = 0.543. Our analysis however indicates stable 
2D rolls for nonzero P, with the range of r p — 1 decreasing 
rapidly with the lowering of P. Earlier Krishnamurti [14] 
had found stable 2D rolls for mercury (P = 0.025) near 
the onset (R sa (2.3 ± 0.1) x 10 3 ) in her experiments with 
no-slip boundary conditions, thus indicating that P c pre- 
dicted by Busse and Bolton [5] is an overestimate. Re- 
garding lu, Busse [4] predicts that lu ~ P in contrast to 
our finding that lu ~ P 2 near the onset of convection. In 
a related work on convection, Mercader et al. [25] have 
shown that the relationship between lu and P depends 
critically on the aspect ratio and they have reported both 
linear and quadratic dependence. These issues require a 
more careful and detailed theoretical investigation. 

A comparison of the above results with the bifurcation 
diagram of zero-P convection of Pal et al. [1] reveal strong 
similarities. The bifurcation scenario from SQ to SQOR 
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Fig. 5: (a) Phase space projection of the chaotic attractor Ch3 
on the Wioi-Won plane obtained from DNS with P = 0.005 
and r = 1.0023; (b) Phase space projection of the chaotic at- 
tractor Chi on the Wioi-VFon plane obtained from DNS with 
P — 0.005 and r = 1.0068. The insets show an enlarged view 
of the boxed region. 

for the larger r values are topologically equivalent. We 
observe two major differences between the bifurcation di- 
agrams of zero-P convection and low-P convection: (a) no 
chaos for P > 0.02 near the onset; (b) existence of stable 
2D rolls and associated ASQ and OASQ patterns near the 
onset. For P = 0.02 we find only stationary and time- 
periodic solutions near the onset which is in agreement 
with the results of Meneguzzi et al. [18] for P = 0.025 and 
Thual [19] for P = 0.02 and 0.2. However, for P < 0.005 
we observe chaotic attractors Chi, Ch2 and ChS between 
OASQ (closer to r = 1) and SQOR. The stationary 2D 
rolls and associated ASQ and OASQ patterns are observed 
in a narrow window between the conduction state and the 
chaotic states. 

In conclusion, we perform a bifurcation analysis for 
low-P convection near the onset using DNS and a re- 
lated low-dimensional model. The results of the low- 
dimensional model are in good agreement with those of 
DNS since the model is derived using DNS by choosing 
most of the relevant modes. We have chosen aspect ratio 
of T x = Ty = 2y/2. We observe that the bifurcation di- 
agram for low-P convection is very similar to the zero-P 
convection reported by Pal et al. [1], except near the on- 
set of convection where 2D stationary rolls, and station- 
ary and oscillatory asymmetric squares are observed for 
nonzero Prandtl numbers. The range of Rayleigh num- 
bers for which 2D rolls and associated ASQ and OASQ 
are observed shrinks rapidly (~ P 2 ) as P is decreased. 
This result is in qualitative agreement with the results of 
Bussc [4]. For P < 0.0002, the range of reduced Rayleigh 
numbers for which the stationary 2D rolls could be ob- 
served is too narrow (< 10~ 7 ) to be observed in experi- 



Fig. 6: Bifurcation digram of the low-dimensional model for 
P = 0.0002 in the range 0.95 < r < 1.247. The color conven- 
tion is same as that for P = 0.005 (Fig. [3]). A zoomed view of 
the bifurcation diagram for the chaotic regime is shown in the 
inset where the x axis is chosen as log(r — 1) to highlight the 
behaviour near r = 1. 

ments or in DNS. Our results predict the critical Prandtl 
number P c to be much lower than that predicted by Busse 
and Bolton [5]. 

In the present letter we also report the origin of chaos 
near the onset of convection for low-P convective flows. 
Here chaos appears through a series of crisis. The same 
set of chaotic attractors also appear in zero-P convec- 
tion [1]. Thus the bifurcation diagrams for low-P convec- 
tion (P < 0.005) have significant similarities with zero-P 
convection, and zero-P convection is an asymptotic limit 
of low-P convection. Our bifurcation analysis provides 
very useful insights into the origin of patterns and chaos 
for low-P convection. An extension of the present work 
to different aspect ratios, and to low-dimensional models 
with a larger set of modes is in progress. 
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